Initialization of NMR Quantum Registers using Long-Lived Singlet States 
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An ensemble of nuclear spin-pairs under certain conditions is known to exhibit singlet state life- 
times much longer than other non-equilibrium states. This property of singlet state can be exploited 
in quantum information processing for efficient initialization of quantum registers. Here we describe 
a general method of initialization and experimentally demonstrate it with two-, three-, and four- 
qubit nuclear spin registers. 

PACS numbers: 



I. INTRODUCTION 

It has been established theoretically that quantum sys- 
tems are far more capable than their classical counter- 
parts for certain computational tasks [1] . Although vari- 
ous types of quantum systems are being explored, the re- 
alization of a general purpose quantum processor remains 
a technological challenge 2] . Proof-of-principle type ex- 
perimental demonstrations have been carried out for sev- 
eral problems including factorization, unsorted database 
search, and quantum simulations using Nuclear Magnetic 
Resonance (NMR) 0, 0| . Nuclear spins in an ensemble 
of molecules present a convenient architecture to simu- 
late a quantum register. The practical demonstrations of 
Quantum Information Processing (QIP) using such reg- 
isters are greatly benefited by the long coherence times of 
nuclear spins and the well developed control techniques 
of NMR [111]. 

In order to carry out information processing, a quan- 
tum register must satisfy a set of criteria laid out by 
DiVincenzo [f|. An important criterion is the ability to 
precisely initialize the register to a desired ket of the com- 
putational basis. Nuclear spins in room temperature and 
at ordinary magnetic fields exist in a highly mixed state 
and therefore preparing a pure state is not straightfor- 
ward [6-8]. The spin temperature can nevertheless be 
reduced by using parahydrogens @ or by using Dynamic 
Nuclear Polarization In future, either or both of 

these techniques may be available for preparing NMR 
quantum registers into almost pure states The ex- 
isting approach for initializing NMR registers is however 
based on specially prepared mixed states known as 'pseu- 
dopure states', which are isomorphic to pure states for 
several computational problems [TJ [l2[ • 

Consider an ensemble of identical molecules each hav- 
ing n spin-1/2 nuclei in a magnetic field. The Zeeman 
Hamiltonian rl z = ^oHi 1S characterized by the fre- 
quency Vq of Larmor precession, and the z-component 
of the spin angular momentum operator P z of spins 
j = 1 • • • n 0. In NMR-QIP, the eigenstates |±l/2) of U{ 
are labeled as |0) and |1) states of a qubit, and the multi- 
spin eigenbasis {|00 • • • 00), |00 • • • 01), ■ • • } is treated as 



the computational basis. In thermal equilibrium at or- 
dinary room temperature T, the Boltzman factor kT is 
much larger than the Zeeman energy gaps, so that the 
density matrix can be expanded as 



Pcq 



= 2-"e"^/ feT « 2-"(l + p A ). 



(1) 



Here identity 1 corresponds to a background of uniformly 
populated levels and the trace-less part p& — Y]. €jl j z is 
known as the deviation density matrix. The dimension- 
less numbers Ej — —hvL/kT have magnitudes ~ 10~ 4 
for protons in currently available magnetic fields at or- 
dinary room temperatures. The deviation density ma- 
trix represents the unequal population distribution (over 
the uniform background) which lead to the observable 
magnetizations. Thus preparing a pure ground state i.e., 
|000 • • • 0) in an NMR quantum register is rather difficult. 
Nevertheless, Cory et al 11 1 and Chuang et al [l2| have 
independently pointed out that often a pure state IV'KV'I 
can be simulated by the pseudopure state 



Ppps 



[(l-e')l + 2V|^|] 



(2) 



Here e' is a measure of the magnetization retained in the 
pseudopure state and it usually gets halved with every 
additional qubit jl3| . The unit background is invariant 
under the Hamiltonian evolution, does not lead to NMR 
signal and is ignored pd| . Thus the equilibrium den- 
sity matrix of a single spin-1/2 nucleus is always in a 
pseudopure state. Initializing a multi-spin system into 
a pseudopure state however is essentially a non unitary 
process Qj. Several methods have earlier been proposed 
and they involved averaging of magnetization modes over 
the sample space (called 'spatial averaging' fill ]) or over 
spin space (called 'logical labeling' [12|, fl5j | ) o r over sev- 
eral transients (called 'temporal averaging' |l6(). In some 
cases subsystem pseudopure states are easier to prepare 
either by transition selective pulses 
selection using pulsed field gradients 



17] or by coherence 
181 ] , but these meth- 



ods invariably result in loss of a qubit for further compu- 
tation. 

In the next section we propose a different approach 
that exploits long life-times of certain special states called 
'singlet states'. The following section details the exper- 
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imental demonstrations on model systems consisting of 
two, three and four-qubit NMR registers. 



detection of singlet states and then elucidate the initial- 
ization of quantum registers. 



II. SINGLET STATES AND REGISTER 
INITIALIZATION 

A. Singlet States 

The Hamiltonian for an ensemble of spin- 1/2 nuclear 
pairs of same isotope, in the RF interaction frame, can 
be expressed as 



H cS = h 



Au _ Au 

2 z 2 z 



JI 1 ■ I 2 + vail 



(3) 



Here the RF frequency is assumed to be at the mean 
of the two Larmor frequencies, and Av, J and i>i2 corre- 
spond respectively to the difference in Larmor frequencies 
(chemical shift difference), the scalar coupling constant 
and the RF amplitude (all in Hz). In the limiting case 
of Av — > 0, the system is said to have magnetic equiva- 
lence, and the singlet state |So) = (|01) — |10))/%/2, and 
the triplet states |T X ) = |00), \T ) = (|01) + |10))/V2, 
and \T-\) = |11) form an orthonormal eigenbasis of the 
internal Hamiltonian H c ^ = hJI 1 ■ I 2 |[. 

The non-equilibrium nuclear spin states decay toward 
equilibrium state ([T]) with a time constant called 'spin- 
lattice' relaxation time constant Ti 7]. Since most of 
the NMR experiments involve preparation and detection 
of non-equilibrium spin states, it was generally accepted 
that the duration of any NMR transient is limited by T\ , 
although the theoretical limit of transverse relaxation is 
2Ti Q. Levitt and co-workers demonstrated that un- 
der the Hamiltonian Ti.^ the decay constant of singlet 
state | So) is much larger than T±, and hence called it a 
long-lived state [ltl [20| . The phenomenon is akin to the 
decoherence-free-subspace (DFS) which is well known in 
QIP [2l| . Detailed theoretical analyses of singlet state de- 
cay have been provided by Levitt and co-workers [22|, HH 
and by Karthik and Bodenhausen (24|. The long life 
times of singlet states have been attributed to the fact 
that the singlet states are immune to the dominant in- 
tramolecular dipole-dipole relaxation mechanism, which 
is symmetric with respect to the exchange of spin states. 

Experimentally, sing let decay constants upto 36Ti 
have been reported [25[, and in another instance singlet 
lived till about half an hour (2||. Overcoming the Ti 
limit has motivated a number of novel applications in- 
cluding studying slow molecular dynamics and transport 
process es [ 25l . |27| , precise measurements of NMR inter- 
actions [28J , the transport and storage of hyper polarized 
nuclear spin ord er |29l - l34j . homogeneous line narrowing 
in spectroscopy |35j|. and determining molecular torsion 
angles [3(3] . Recently we had reported that high-fidelity 
singlet states can be easily prepared and characterized in 
two-spin systems under various experimental conditions 
[37| . In the following we describe the preparation and 



B. Preparation and detection of singlet states 

The state |So)(So| — |T )(T | can be easily prepared 
from the equilibrium density matrix l\ + 1 2 by using the 
propagator 



Ul* = e -iWl-il) e -^l 



■(4) 



With certain approximations the above propagator can 
be constructed from the Hamiltonians similar to the one 
in © [38]. The equivalence Hamiltonian may be 
realized by suppressing the chemical shift Av either by 
vanishing the Zeeman field [l9| . or by using a 'spin- lock' 
[20l l39j . In this work we employ the latter technique. 
The spin-lock may be achieved by applying a long low- 
power unmodulated RF [2(J, or by specially designed 
phase modulated sequence which were originally used for 
spin-decoupling [39j |. 

The singlet states by themselves are inaccessible to 
macroscopic observables, but can be indirectly detected 
by removing the equivalence and transforming to ob- 
servable single quantum coherence using the propagator 

MM 



u 



(5) 



Alternatively, a more detailed and quantitative analysis 
of singlet states may be carried out using density matrix 
tomography (37j . 



C. Initializing NMR Registers 

The circuit for initializing a 2-qubit NMR register via 
singlet states is shown in FigQJi. An initially imperfect 
singlet density matrix gets purified during the spin-lock 
period as a result of the long life time, while the artifact 
coherences are destroyed by relaxation process as well as 
the inhomogeneities in the spin-lock itself. There exist 
optimal spin-lock conditions at which one obtains singlet 
states with high fidelity [37[. Once the singlet state is 
prepared with high fidelity, the conversion |So) — > 1 01) 
can be easily achieved by the propagator 



e 2 * 



(0) 



Finally a pulsed field gradient G z can be used to destroy 
the residual single and multiple quantum coherences. If 
necessary, other pseudopure states can be obtained sim- 
ply by applying NOT gates. 



This procedure of initialization can be extended to 
multiqubit systems, since it is known that a spin-pair 
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spin-lock of multiple singlet pairs can be achieved using 
sophisticated modulated RF sequences as described in 
the next section. 
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FIG. 1: The circuit diagrams for initializing (a) 2-qubit, (b) 
3-qubit, and (c) n-qubit registers. In (c), the cNOT gates 
with open circles correspond to NOT operation if the control 
is and identity if the control is 1. The h-gate corresponds 
to pseudo-Hadamard: |0) A (|0) - |1»/V2 and |1) b 
|X»/V2. 



(|0> 



may exhibit long-lived singlet state if it is situated suf- 
ficiently away from other spins [36[ . The procedure of a 
3-qubit initialization is shown in FigQJ>. We first prepare 
a two-qubit singlet and apply a cNOT gate on qubit-2 
(controlled by qubit-3). Subsequent spin-lock and 
gate on qubits 1 and 2 initializes a three qubit system 
into |010) state. The circuit can be understood as fol- 
lows. After preparing the singlet on qubits 1 and 2, the 
third qubit remains in a mixed state with a probability 
Po of being in state |0) and a probability pi of being in 
state |1). The cNOT gate transforms the mixed state 
according to: 



cNOT 



|5o)(5 |®(po|0)(0|+ Pl |l)(l|) 
PolSbXSbl «» |0><01 +jpi|^-><0-| O |1><1| 



(7) 



First term in the output is singlet on qubits 1 and 2 
and therefore survives during the second spin-lock where 
as the other term consisting of the Bell state \<t>~) — 
(| 00) — |ll))/\/2 decays fast. The singlet is ultimately 
transformed into |01) by U\ ,rL and thus we obtain 1 010) 
pscudopure state with a good approximation. 

A general scheme for the initialization of an n-qubit 
NMR quantum register is described in FigQJi. The circuit 
can be analyzed in a similar way as in the 3-qubit case. 
The fact that only nearest-neighbor interactions are used 
is highly advantageous in practice. Experimentally, the 



III. EXPERIMENTS 

The following experiments are carried out in Bruker 
500 MHz NMR spectrometer at an ambient tempera- 
ture of 300 K. We used strongly modulated pulses for 
designing high fidelity local gates as well as cNOT gates 
[HSU. The spin-lock was achieved by WALTZ-16 - a 
phase modulated RF sequence, which is routinely used 
in broadband spin decoupling [37|. Spectra correspond- 
ing to pseudopure states are obtained by linear detection 
scheme using small flip angle RF pulses. Since the pseu- 
dopure states have one energy level more populated than 
the equal distribution in all others, the spectrum should 
consist ideally of only one transition per qubit in each 
case. Quantitative analyses of the pseudopure states are 
carried out using extended versions of density matrix to- 
mography described in reference [13] . Finally, the success 
of the experimental state p in achieving a target pseudop- 
ure state /Opp S is measured by calculating the correlation 
401. 



(P 



trace [p ■ p pps ] 



pps/ 



trace [p 2 ] ■ trace [pp ps ] 



(8) 



Often only the diagonal elements of the density matrices 
are relevant and in such cases, the 'diagonal correlation' 
can be expressed by replacing all the operators in the 
above expression by their diagonal parts. In the following 
we describe the individual cases of two-, three- and four- 
qubit registers. 



A. 2-qubit register 

The two-qubit system, Hamiltonian parameters, and 
the corresponding pseudopure and the reference spectra 
are shown in Figj2ja-d). As shown in FigJT}i, the exper- 
iment involved preparing singlet using U 1 ' , followed by 
RF spin-lock with amplitude 2 kHz and duration 12.4 s, 
which are optimized for high singlet content [371 ] . The 
decay constant for singlet state was 16.2 s approximately 
three times the T\ values of the two spins. The singlet 
is then converted into 1 01} pseudopure state using f7 2 ' . 
A final gradient pulse served to destroy the artifact co- 
herences. The bar plots showing the real and imaginary 
parts of the theoretical and experimental density matrix 
are shown in Figl^e-h). A very high correlation of 0.995 
is obtained with 1 01} pseudopure state. 



For some quantum algorithms a good starting point 
may be the singlet state itself, which can be extracted 



4 



(a) ^ m (b) 



(c) 



(d) 



100 50 




1 


2 




96 


4 . 


1 




-96 


2 



-50 -100 Hz 




directly after the spin-lock shown in FigQJi. The real part 
of the experimental singlet density matrix shown in FigJ5J 
has a correlation of 0.991. If necessary, initialization to 
other Bell states, starting from the singlet state, can be 
carried out easily: 



|So) 
|So> 
I So) 



h/> + ) = (|01> + |10»A/2, 
= (|00) + |11))/V2, 
|0-> = (|OO)-|11»/V2. 



(9) 



The z-rotation in the above propagators can be imple- 
mented by simply chemical shift evolution for a period 
1/(2A^), and the qubit selective x-rotation can be im- 
plemented by using a strongly modulated pulse. The ex- 
perimental density matrices corresponding to these Bell- 
states have respective correlations 0.987 (FigJ2J), 0.982 
(Figl2]0, and °- 968 (FigH). 



B. 3-qubit register 




FIG. 2: The molecular structure (a) and Hamiltonian pa- 
rameters (b) of 5 -bromothiophene-2-carbaldehyde dissolved 
in CD3OD, forming a homonuclear two-qubit register. In (b) 
diagonal and off-diagonal elements correspond to the chem- 
ical shifts and the scalar coupling constant respectively (in 
Hz). The 1 H spectra correspond to the pseudopure state (c) 
and the equilibrium mixed state (d). The barplots (e-h) corre- 
spond to the real (e,g) and imaginary (f,h) parts of theoretical 
(e,f) and experimental (g,h) pseudopure 1 01) state. The re- 
maining barplots correspond to the experimental real parts of 
\So) (i), |^ + ) G), 10+) (k), and |«T> (!)• 



The three-qubit system, Hamiltonian parameters and 
the corresponding pseudopure and the reference spectra 
are shown in FigJ3ta-d). The decay constant for singlet 
state of spins 1 and 2 was about 18 s, approximately three 
times of their T\ values. The pseudopure state was pre- 
pared using the circuit shown in Fig0J>. Each of the two 
spin- locks were achieved using 6.3 s of 500 Hz WALTZ- 
16 modulations. The cNOT gate was implemented using 
a 14 segment strongly modulated RF pulse of duration 
approximately 60 ms and of fidelity 0.96. The bar plots 
showing the real and imaginary parts of the theoretical 
and experimental density matrix are shown in FigJ3](e-h). 
The correlation of the experimental density matrix with 
the theoretical pseudopure state is |010) is 0.952. The 
correlation is smaller compared to the two-qubit case, 
mainly due to the errors in the cNOT gate. Neverthe- 
less, the diagonal correlation is as high as 0.983. 



C. 4-qubit register 

The four-qubit system, Hamiltonian parameters and 
the corresponding pseudopure and the reference spectra 
are shown in Fig0] The pseudopure state was prepared 
using the circuit shown in FigfTJ;. The singlet decay con- 
stants were about 6 s, approximately twice the Ti values 
of the individual spins. We were able to carry out si- 
multaneous spin-lock of two singlet pairs and initialize a 
four-qubit register. The two spin-locks were achieved by 
2 kHz WALTZ-16 modulations of durations 2 s and 4.5 s 
each. The two cNOT gates were made of 20 segments, ap- 
proximately 61 ms duration and of fidelities about 0.94. 
The 10 segment h-gate was about 8.2 s long and of fidelity 



5 



(a) 



(b) 



1 


2 


3 




38 


. 9 


17 . 9 


1 




-38 


11 . 7 


2 






-244 


3 




50 40 Hz -30 -40 Hz -230 -240 -250 Hz 




110 



(e) 11 °^ooo 010 




110 



(f ) 110 ^„„010 



000 





000 
010 
100 

(h) 110 



110 



000 



,010 



FIG. 3: The molecular structure (a) and Hamiltonian pa- 
rameters (b) of acrylonitrile dissolved in CDCI3, forming a 
homonuclear 3-qubit register. In (b) diagonal elements are 
chemical shifts (in Hz) and the off-diagonal elements are the 
scalar coupling constants (in Hz). The 1 H spectra correspond 
to the pseudopure state (c) and the equilibrium mixed state 
(d). The bar plots are showing the real (e,g) and imaginary 
(f,h) parts of theoretical (e,f) and experimental (g,h) pseu- 
dopure |010)(010| state. 



0.98. Complete tomography of a 4-qubit density matrix is 
a laborious task. After the preparation of the pseudopure 
state, the non-zero quantum off-diagonal elements are ef- 
ficiently destroyed by the final gradient pulse. Since only 
the diagonal elements are of main interest, we have car- 
ried out the four-qubit diagonal tomography 41]. The 
bar plot showing the diagonal part of the experimental 
density matrix is shown in FigHJ:. The diagonal corre- 
lation is estimated to be approximately 0.97 ± 0.01 with 
1 1001) pseudopure state. The first pair has collapsed to 
1 10) state instead of |01), due to an additional 180 degree 
pulse that was applied on qubits 1 and 2 for refocusing 



purposes during U 2 
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FIG. 4: The molecular structure (a) and Hamiltonian param- 
eters (b) of aspirin dissolved in CD3OD, forming a homonu- 
clear 4-qubit register. In (b) diagonal elements are chemical 
shifts (in Hz) and the off-diagonal elements are the scalar cou- 
pling constants (in Hz). The barplot in (c) displays the diag- 
onal elements of the density matrix obtained by tomography 
of the pseudopure 1 1001) state. The 1 H spectra correspond 
to the pseudopure state (d) and the equilibrium mixed state 
(e). 



IV. CONCLUSIONS 



We have demonstrated that robust initialization of 
NMR quantum registers are possible via long-lived sin- 
glet states. It is hard to initialize proton-based NMR 
registers using traditional methods. As a result many 
popular NMR registers were based on carbon spins or a 
combination of protons and carbons which permitted ini- 
tialization by traditional methods. The proposed method 
of using long-lived singlet states enables us to initialize 
larger registers even though the long range interactions 
are weak. Molecules are of interest where in the inter pair 
dipolar couplings are sufficiently weak to keep the sin- 
glet states long-lived, while the covalant bond mediated 
scalar interactions among the nearest neighbor spins are 
sufficiently strong. Since para hydrogens naturally ex- 
ist in singlet states, the method can be applied directly 
to the initialization of registers based on parahydrogens. 
While the register initialization is the first application of 
long-lived singlet states for QIP, more applications, like 
for example enhancing the memory of registers, may be 
realized in future. Similar techniques may be used for 
multi-qubit initialization in non-NMR systems exhibit- 
ing long-lived singlet states. 
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